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Abstract 

The dynamics of "dipolar particles", i.e. particles endowed with a four- vector mass dipole 
moment, is investigated using an action principle in general relativity. The action is a specific 
functional of the particle's world line, and of the dipole moment vector, considered as a dynamical 
variable. The first part of the action is inspired by that of a particle with spin moving on an 
arbitrary gravitational background. The second part is intended to describe, at some effective 
level, the internal non-gravitational force linking together the "microscopic" constituents of the 
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dipole. We find that some solutions of the equations of motion and evolution of the dipolar 
particles correspond to an equilibrium state for the dipole moment in a gravitational field. Under 
some hypothesis we show that a fluid of dipolar particles, supposed to constitute the dark matter, 
reproduces the modified Newtonian dynamics (MOND) in the non relativistic limit. We recover the 
main characteristics of a recently proposed quasi-Newtonian model of "gravitational polarization" . 
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I. INTRODUCTION 



A. Astrophysical motivation 

It is recognized (see Refs. [3, 0,(1, 0] for reviews) that the modified Newtonian dynamics 
or MOND, which has been proposed by Milgrom [H, @, @] as an alternative to dark matter, 
works extremely well at predicting the form of rotation curves of galaxies from the observed 
distribution of stars and gas. In addition MOND explains naturally the relation between the 
luminosity of galaxies and their asymptotic rotation velocity — the so-called Tully-Fisher 
law [8]. Though there might be some examples where MOND does not fully account for 
the observed kinematics of galaxies 0], overall the fit achieved by MOND of the rotation 
curves of most galaxies is very impressive and calls for a possible physical explanation. On 
the other hand we know that the mass discrepancy of clusters of galaxies is not completely 
accounted for by MONDjlOj, and that at the cluster scale there is still empirical evidence 
for unseen dark matter 

In the usual interpretation, MOND is viewed as a modification of the fundamental law of 
gravity, without the need of dark matter. Several relativistic extensions of MOND, sharing 
this view of modifying the sector of gravity, postulate the existence of supplementary fields 
associated with the gravitational force, in addition to the metric tensor of general relativ- 
ity T2I, 13, 14, 15, ly]. So far the most developped of the relativistic MOND theories is the 
tensor-scalar- vector theory of Bekenstein & Sanders [3, 15, 16 . 

At this stage we seem to face two alternatives to the issue of dark matter: 

(1) Either accept the existence of cold dark matter particles, e.g. predicted by super- 
symmetric extensions of the standard model of particle physics (see [l7[ for a review). 
However these particles are yet to be discovered, and the simplest models of cold dark 
matter fail to reproduce in a natural way the flat rotation curves of galaxies ^ . 



(2) Or postulate an alteration of our fundamental theory of gravity (namely MOND and 
its relativistic extensions). But the motivation for altering the law of gravity is ad 
hoc, and one can argue that the relativistic MOND theory (l4l . flR 16| does not in 



fact explain the phenomenology of MOND — at least untill the extra fields find a 
well-motivated explanation coming from fundamental physics. 

In the present paper, following [l8[ 1 , we propose a third alternative. Namely we: 

(3) Keep the standard law of gravity, i.e. general relativity and its Newtonian limit, but 
we add to the distribution of ordinary matter some specific form of dark matter in 
such a way as to naturally explain MOND. The dark matter consists of "polarization 
masses" associated with a medium of dipole moments aligned in the gravitational field 
of ordinary matter. (But, in this paper, we shall leave aside the problem of clusters of 



galaxies [10|, 111].) 



Our basic motivation is that MOND can be naturally interpreted as resulting from an 
effect of gravitational polarization of a dipolar medium. Paper I argued, on the basis of 
a simple quasi Newtonian model, that the polarization tends to enhance the magnitude 
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18| will be referred to as paper I. 
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of the gravitational field of ordinary galaxies in a way consistent with MOND. This effect 
constitutes the gravitational analogue of the electric polarization of a dielectric material 
(whose atoms can be modelled by electric dipoles) by an applied electric field [l9(. Thus 
the phenomenology of MOND results from the non-standard influence of the dark matter 
on ordinary matter. 

Arguably the proposal (3) does not really provide an explanation for MOND. Indeed the 
models of paper I and this paper are only effective — the fundamental nature of the dipolar 
particles (i.e. their "internal" structure) will not be elucidated. However the proposal is 
conceptually simple and fits naturally with the MOND phenomenology of the flat rotation 
curves and the Tully-Fisher empirical relation. In addition the model cannot be said to be 
ad hoc because it invokes the (gravitational analogue of the) well-known physical mechanism 
of polarization by an external field. 

The quasi-Newtonian model of paper I is based on a microscopic description of the dipole 
moment using negative gravitational- type masses (or gravitational charges). As a result 
the motion of dipolar particles in this model violates the equivalence principle. In the 
present paper we elaborate a model of dipolar particles and gravitational polarization in 
the standard general relativity theory, 2 without negative (passive) gravitational masses and 
consistent with the equivalence principle. Consequently, we shall find that the equations 
of motion of the dipolar particles in the non-relativistic (NR) limit of the present model 
are different from those of paper I. However we shall recover the main characteristics of 
the quasi-Newtonian model of paper I, and notably the interesting connection with the 
phenomenology of MOND. 



B. Concept of dipole moment in general relativity 

In theories satisfying the equivalence principle the mass dipole moment of a mass distri- 
bution, say 7?* = ^m g x l (adopting a Newtonian picture), is proportional, by equivalence 
between the gravitational and inertial masses, m g = m ; , to the position of the center of mass 
of the system, C l = ^rriiX 1 /M- x (where Mj = Yl m i)- Similarly, the current dipole moment, 
/ig = m g £lJ , is equivalent to the spin angular momentum S l = X] m i £ &v k . Defin- 
ing the mass dipole moment of a particle (supposed to be composed of some sub-particles) 
is a priori delicate because one faces the problem that 7r g = in the center-of-mass frame 
where C % = by definition; however the current dipole moment /x g equivalent to the spin 
S 1 is admissible. Thus, while the notion of a particle carrying a mass dipole moment seems 
to be possible only at the price of violating the equivalence principle (like in the model of 
paper I), the notion of a particle endowed with a current dipole moment or spin is perfectly 
legitimate. 

Particles with spins have been the subject of many fundamental works in general relativ- 
ity 21, 22, [23|, 0, 25, |26|. with app lications to the problem of dynamics of spinning black 
holes in binary systems [27l. [28l . |29| . [30I . 3l| . The general formalism has been encapsulated 



in an action principle for a particle with s pin moving on an arbitrary gravitational back- 
ground. The action, due to Bailey & Israel 26[ , is derived by expressing the standard action 



for a non-spinning particle (i.e. the integral of the proper time), as defined with respect to 



2 The concept of gravitational polarization at the quadrupolar order in relativistic gravity theories has been 
investigated in Ref. Q. 
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some "eccentric" world-line located at the edge of some composite particle, in terms of the 
quantities belonging to a different world-line, which is viewed as the "reference" or "central" 
world-line located at the center of the particle. Crucial to the formalism is the vector separa- 
tion between the two world-lines, defined as the gradient of the geodesic distance separating 
them. 3 The spin of the particle is then given by the anti-symmetric product between the 
latter separation vector and the particle's linear momentum. Hence the separation vector 
is the "lever arm" associated with the spin angular momentum, as defined with respect to 
the fiducial central world-line. The point for our purpose is that such lever arm appears 
essentially to be the analogue of a mass-type dipole moment. We thus see the possibility 
of describing dipolar particles in general relativity by starting from (some variant of) the 
Bailey-Israel action [26l |. 

The previous paradox concerning the dipole moment which should vanish in the frame of 
the center-of-mass will be solved because (roughly speaking) there are two notions of dipole 
moments. The first moment is the one we mentionned earlier, 7r g , but let us call it now the 
"passive" dipole moment 7r p to emphasize the fact that the gravitational mass it contains 
is really the passive gravitational mass m p = m g (i.e. that mass which enters in the right- 
hand-side of the law of motion, in factor of the gravitational field). Assuming as before that 
the dipolar particle is composed of sub-particles, its passive dipole moment 7r p ~ Yl m p x % 
will indeed be zero in the particle's center of mass, by equivalence between the inertial and 
passive gravitational masses of the sub-particles, m p = m\. However the dipole moment 
we shall consider in this paper is different: this is the moment parametrizing the dipolar 
part of the stress-energy tensor T^ v of the particles (notably its 00 component or energy 
density). In general relativity the stress-energy tensor is the source for the gravitational 
field, so we can rightly refer to this moment as the active dipole moment. In conclusion, we 
shall have 7r p = by the equivalence principle, but we still have at our disposal the "active 
gravitational" version of the moment, call it 7r*, which enters the stress-energy tensor of 
the particles, say through a term of the form ~ —di(ir\8) in the energy density at the NR 
approximation (where 5 is the Dirac function). In the following we shall investigate a specific 
model for the relativistic dynamics and evolution of particles endowed with a four-vector 
dipolar moment 7r M = n£ of the active type. 4 

By analogy with the model of paper I we shall denote by 2m the total inertial mass of 
the dipolar particle, equivalent to its total passive gravitational mass, i.e. M; = M p = 2m. 
However, this equivalence does not mean that the motion of the dipolar particle is geodesic. 
Indeed, motivated by paper I, we shall introduce a force F^, which will be considered as 
"internal" to the dipolar particle, and is aimed at "stabilizing" the dipole moment imbedded 
in an exterior gravitational field. The force F 11 has a non gravitational origin, and will 
derive from a scalar potential function V in the action. Because of the presence of this 
force reflecting its internal structure, the dipolar particle is not a test particle: its motion is 
not geodesic, its four-acceleration is non zero; the particle can be thought of as a "rocket" , 
self- accelerated by the internal force F^. We shall find an approximate solution in which 
F^ accelerates the particle in such a way as to compensate for the local gravitational field, 



3 The geodesic distance is nothing but the world function in Synge's formalism [32[ ; it is a bi-scalar in the 
general theory of bi-tensors . 

4 Hereafter we no longer mention that this moment can be regarded as the active one. Note that the names 
we are giving here and below for different types of masses and dipoles are useful for the discussion but 
are not required in the present formalism, which is entirely based on the action (|2.1|) . 
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so that the dipolar particle stays essentially at rest in a gravitational field (like in the model 
of paper I). 

The monopolar part of the stress-energy tensor defines what can be regarded as the 
particle's active gravitational mass M a . Intuition from paper I would lead to expect that for 
a dipolar particle M a is zero, i.e. the particle does not generate any monopolar gravitational 
field. In the present model we shall find that this is not possible: M a cannot be zero, and 
in fact we shall have M a = V/c 2 , where V is the potential function from which derives the 
internal force F^. However, thanks to the explicit factor 1/c 2 it contains, M a turns out to 
be very small since it vanishes in the NR limit, i.e. M a = 0(c~ 2 ) when c — ► oo. Thus, the 
dipolar particle will indeed be "purely dipolar" in the NR limit. 5 

Applying the standard general relativistic coupling to gravity, we prove that the Einstein 
field equations reduce to a Poisson equation with a dipolar source term in addition to the 
density of ordinary matter, having the structure ~ — di(ir)_5) in the limit c — > oo (where ir\_ 
refers to an appropriate orthogonal projection). We are then close to the MOND equation; 
to recover MOND it suffices that the dipole moment be aligned with the gravitational field, 
and polarized in a certain way. In the present paper we shall find, under some hypothesis, 
some solutions which correspond to gravitational polarization, i.e. we shall show that the 
polarization scenario yielding MOND is consistent with our equations. There will be essen- 
tially a one-to-one correspondence between the internal potential scalar function V entering 
the action and the Milgrom 0, @, 0] function /x (which is linked in the interpretation of 
paper I with the gravitational susceptibility x °f the dipolar medium by \i = 1 + x). 

The paper is organized as follows. We first deal with the relativistic model based on the 
action presented in Sect. Ill Al The general equations of motion and evolution, and the stress- 
energy tensor, are derived in Sect. Ill Bt In Sect. Ill CI we restrict ourselves to a particular 
solution which corresponds to some "equilibrium" state for the dipole moment. For that 
solution we obtain in Sect. MI Al the dynamics of dipolar particles in the NR approximation 
where c — ► oo. By neglecting the tidal gravitational field we find in Sect. MI Bl that when 
the fluid of dipole moments is at rest with respect to the ordinary matter distribution the 
medium is polarized and the MOND equation follows. Finally, we investigate in Sect. IIVI 
(still in the NR limit) the case of a stationary fluid of dipole moments aligned in the central 
gravitational field of a point mass. Section |V] summarizes and concludes the paper. In the 
Appendix |X] we show on general grounds how to vary a dipolar action functional. 



II. RELATIVISTIC DYNAMICS OF DIPOLAR PARTICLES 



A. Action principle 



The starting point of our approach is the formalism for spinning particles in general 
relativity (2ll . 22, 23, 24, 25, 26] and especially the Bailey-Israel action [26]. In this formal- 
ism, the equations of motion of the spinning particle (moving on an arbitrary gravitational 
background) are derived from the extremum of the action with respect to an infinitesi- 



5 Notice that the fact that the particle carries some different active and passive gravitational masses, M a 
and Mp, does not contradict the principle of conservation of the linear momentum (and the law of action 
and reaction). Indeed the particle's stress-energy tensor contains in addition to the monopolar part a 
dipolar contribution (parametrized by 7r M ) which ensures the satisfaction of the usual conservation laws. 
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mal displacement of the particle's world-line; however the spin variable itself is not varied. 
Instead, one postulates some condition for the evolution of the spin, called the spin supple- 
mentary condition (SSC), which essentially tells that the spin is parallel transported along 



the particle's world-line (see e.g. [30|, [3l|]). In the present paper we shall interpret the "lever 
arm" associated with the spin angular momentum, in the formalism of spinning particles, 
as the mass-type dipole moment carried by the particle. However we shall not impose some 
analogue of spin supplementary condition. Instead we shall promote the dipole moment as 
a dynamical variable, which will have to be varied independently of the variation of the 
particle's position. 

Let the particle follow a time-like world line x M in a four dimensional manifold, and being 
endowed with the dipole moment four vector tt^ (having the dimension of a mass times a 
length). 6 We propose that the particle's dynamics in a prescribed gravitational field g^ v be 
derived from the following matter action, consisting of three terms, 



S=Y* f dr -cJg„ v (2m ^ - tW) (2m u v - ix u ) + ~v(—\ . (2.1) 

^-^J-oo v 4m Vm 



The sum goes over all the particles of this type. We denote the particle's proper time by 
dr = [— g^dx^dx 1 * / c 2 } 1 / 2 , and its four-velocity by w M = dx^/dr (which is normalized to 
g fMU u fJl u u = — c 2 ). The dynamical variables are the space-time position x M (r) and the dipole 
moment 7r M (r), both of them are function of the proper time. 

The first term in Eq. (12. ip is inspired by the known action for the dynamics of spinning 
particles (26J, and describes the inertial properties of the dipolar particle, i.e. its equations 
of motion in the gravitational field. The second term is a kinetic-type term for the dynamical 
evolution of the dipole moment itself, and serves to tell how this evolution will differ from 
parallel transport. The third term of (12.11) is made of a scalar potential function V, supposed 
to describe, at some effective "macroscopic" level, a non-gravitational force that is internal 
to the dipolar particle. 

We denote by a dot the covariant derivative of the dipole moment with respect to the 
proper time, namely 

*> = %r = % + (2.2) 

More generally the dot will always refer to the covariant time derivative D/dr. The potential 
function V in the third term of the action depends on the norm tt± of the orthogonal 
projection of the dipole moment perpendicular to the four velocity, defined by 



7T 



/' 



(2.3) 



where _L{|= + u^u v jc 2 is the corresponding projection operator. The vector (I2.3P is 
space-like and its norm is given by 

vr± = [g^irl] 1/2 = . (2.4) 

The fact that V depends on the orthogonal projection rather than on 7r M itself is crucial 
for the present formalism. 



6 Greek indices take the space-time values /i, v — 0, 1, 2, 3; Latin ones range on spatial values i, j = 1,2, 3; 
the convention for the Riemann curvature tensor is R^ vpa = dpT^ — d a T^ p + ■ ■ ■ ; the Lorentzian signature 
is — h ++; all factors in the speed of light c and Newton's constant G are indicated throughout. 
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The mass m parametrizing the action (12. ip is a certain mass parameter associated with 
the dipole moment, and is such that x = ir±/m represents the typical size of the dipolar 
particle. The potential function V will typically be [such as in Eq. (13.271) below] a quadratic 
function of the particle's size x = 7tj_/m. The mass m can be viewed as the relativistic 
analogue of the mass in the model of paper I (we have already discussed in the introduction 
some heuristic interpretation of 2m). 

An important feature of the action (12.11) is that its first term is given by the norm of a 
space-like vector. It will be useful to keep in mind this vector, so we introduce a special 
notation for it: 

p» = , (2.5a) 



where A = . /_i _ ^ZL + I^L_ . (2.5b) 

V mc z Am^c^ 

This space-like vector, satisfying p 2 = +4m 2 c 2 , will not represent the linear four- momentum 
of the dipolar particle (otherwise the particle would be a tachyon). 7 As we shall see, since 
there are other terms in the action, the linear momentum of the particle will differ from 
and be normally time-like. 

Because of the space-like character of p M we note that the action (12.11) makes no sense 
in the limit where the dipole moment vanishes, since this limit corresponds to a particle 
with an imaginary mass or tachyon. In this sense the dipolar particle described by (12. ip 
exists only through the existence of the dipole moment, contrary to a particle carrying a 
spin, which is described both by a spin and a mass, and reduces to an ordinary point mass 
in the absence of the spin. [We are speaking of classical particles, whose spin is given by 
the classical notion of angular momentum.] Here the mass parameter m is not independent 
from the dipole moment; it is such that it \ = mx where x is the particle's size. 

Unlike the case of a spinning particle [26( , no explicit coupling to the Riemann curvature 
tensor has been introduced at the level of the action (12. ip . Recall that the Riemann tensor 
in the Bailey-Israel action |26| comes from the particular way this action is constructed, by 
conveying the point particle action from some eccentric world-line of some composite particle 
to the physical central world-line. As a result, if one varies this action with respect to the 
spin, one obtains the same equation as with the variation with respect to the position — both 
equations containing the same coupling to the Riemann curvature (at least at linear order 
in the spin). In the present approach, by contrast, the dipole moment will be considered as 
an independent dynamical variable, and we shall demand that the variation with respect to 
the dipole moment provides an independent evolution equation. 

The introduction of the scalar potential function V is motivated by the quasi Newtonian 
model of paper I, consisting of a pair of sub-particles, with opposite gravitational masses ±m 
and positive inertial masses +m (in analogy with the electric dipole in electrostatics), and 
interacting via some non-gravitational force making possible the existence of a stable dipolar 
configuration. [Indeed, free gravitational masses with opposite sign will accelerate apart from 
each other, and consequently must be bound by an internal force able to counteract their 
gravitational repulsion.] In the present relativistic model we shall also need to invoke such 



7 This is in contrast with the Bailey-Israel action for a spinning particle 26], which is given by the norm 
of the particle's linear momentum — a time-like vector as usual. 
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a force — supposed to describe some non-gravitational attractive interaction between the 
constituents of the dipole. This force will derive from the potential V in the action. 

Let us however stress that we shall not use any model for the elementary constituents 
of the dipole, nor shall we invoke the notion of a negative mass. The formalism follows 
completely from the action (12. ip . which contains simply the moment 7r M and a positive mass 
parameter m. As already commented in the introduction, the mass 2m represents the total 
inertial mass M ; of the dipolar particle (in agreement with paper I), but here, contrary 
to paper I, this mass can also be interpreted as the passive gravitational mass M p , i.e. 
M{ = M p = 2m. Later we shall argue that one can also define an active gravitational mass 
M a by the monopole part of the particle's stress-energy tensor, and that this mass turns out 
to be negligible in the NR limit, M a = C(c" 2 ). 8 



B. Equations of motion and evolution 

By the principle of stationary action we require that S be unchanged under the infinites- 
imal variation 5x M (r) and 5tc^{t) in the dynamical variables, with the condition that they 
vanish on the boundary of the region of integration, i.e. &r M (±oo) = and 57r M (±oo) = 0. 
We provide in Appendix |A] below a summary of the way we vary the dipolar action. 

Varying first S with respect to the change in the dipole moment, holding the particle's 
world-line fixed, 5x^ = 0, we obtain an equation taking the form of the force law 

P^ = -2F^, (2.6) 

where we remind that the dot means the covariant derivative with respect to proper time: 
P M = PP M /dr. In the left-hand-side (LHS) we have defined what will turn out to be the 
linear momentum P M of the dipolar particle; it is given in terms of the space-like vector 
defined in Eq. (12.51) by 

p» = p » + tf> . (2.7) 

As we see P M differs from and we shall check later that P M is time-like. In the right- 
hand-side (RHS) of (|2.6|) . P M is the quadri- force derived from the potential V present in the 
action, and is given by 

F» = <^m. (2.8) 
n±_ ax \m / 

Here dV/dx is the derivative of V with respect to its natural argument x = -K±/m. This 
force is proportional to the orthogonal projection of the dipole moment (12.31) . thus it satisfies 
the constraint 

u^ = 0. (2.9) 

We next perform a variation 5x^ of the particle's position, holding the components of the 
dipole moment 7r M "constant" during the displacement of the world-line. For instance one 
can think of the dipole moment vector as being transported parallely along the displacement 
vector Sx^. We find that the variation of the covariant time derivative of the dipole moment, 



Notice the difference with the quasi-Newtonian "microscopic" model of paper I, in which we have M; = 2m 
and M p = M a = 0. This represents the fundamental difference between the model of paper I and the (NR 
limit of the) present relativistic model. 



s 



7r M = Dn^/dr, yields a Riemann curvature term, which can be understood as coming from 
the non commutation of the variational derivative S with the covariant derivative D/dr. 
During the variation we must take into account the contribution due to the orthogonal 
projection operator _L£ present in the term V(ir±/m). All computations done we end up 
with the equation (where Cl^ = DQ^/dr) 

w = -^w pv y^p° . (2.10) 

See Appendix El for a general derivation of this equation. The RHS represents the analogue 
of the famous Papapetrou coupling to the Riemann curvature tensor in the equation of 
motion of a particle with spin [2l|]. 9 The "linear momentum" f2 M in the LHS of (12.101) is 
found to be different from P p and to be given by 

W = ui i -jf, (2.11) 

in which we used again the convenient definition of p M , Eq. (12.51) . and where uo p represents 
another intermediate quantity given by 



u 



I'- 



= — [-?— + V\ --^F". (2.12) 

cr \ Am J mc z 

At this point we observe that the complete dynamics of the dipolar particle is encoded into 
the two equations (12.61) and (I2.10j) . 

Next we obtain the particle's stress-energy tensor by varying the action (12.11) with respect 
to an infinitesimal change in the background metric, Sg^, vanishing at the edges of the 
space-time manifold, when — > oo. Obviously we must take into account all metric 
contributions, including crucially those arising from the Christoffel symbols in the covariant 
time derivative n^, and those coming from -L fJiU = g^ + u^u v /c 2 . The conserved number 
density n of the dipolar particles satisfies the covariant continuity equation 

V v (n/)=0. (2.13) 

By straightforward calculations — see Appendix [A] — we find that the stress-energy tensor 
T pu (with the dimension of an energy density) of the dipolar particles can be expressed in 
terms of the two basic linear momenta and P M as 

T" v = - 7^V p (n tt p P^ - P p n^ u u) ^j . (2.14) 

We readily verify that the covariant conservation law 

V u T^ = 0, (2.15) 

holds as a consequence of the equations of motion (12. 6p and (I2.10p . 

Some physical interpretation follows from the expression of T^ u . It is clear that the 
first term in (12. 14|) takes the form of a monopolar contribution, appropriate for a point-like 



9 In the case of a spinning particle the antisymmetric spin tensor is given in terms of the dipole moment 
variable and the particle's linear momentum by = iMP 1 ^ jm. 



9 



particle having velocity and linear momentum VL^ . We thus see that the particles we 
are considering are not purely dipolar, since they involve also a monopolar contribution 
in this sense. The monopolar piece in the stress-energy tensor will generate a monopolar 
gravitational field via the Einstein field equations. So the mass associated with the linear 
momentum f2 M can be naturally interpreted — since it represents the point-like source of a 
monopolar gravitational field — as the active gravitational mass of the particle [see (I2.3ip 
below for the computation of this mass within a particular solution of the equations] . Simi- 
larly fi M can be referred to as the "active linear momentum" , while P M which enters the law 
of motion (12.61) should rather be regarded as the "inertial linear momentum". [Technically 
f2 M represents the conjugate momentum of and P M the conjugate momentum of 7r M ; see 
Eqs. ( lAlip in the Appendix.] 

The second term in Eq. (12.141) is clearly dipolar, and represents the relativistic gener- 
alization of the quasi-Newtonian density of polarization p po iar = — SjlP in paper I — i.e. 
minus the divergence of the polarization vector II 1 . To emphasize the dipolar character of 
this term we rewrite the stress-energy tensor as 

= n flK v) - c 2 V p fP^ , (2.16) 

where U pfiu can be called the polarization tensor and is given by 



n 



2mc 2 



7T 



u v) . (2.17) 



The gravitational field generated by the distribution of dipolar particles will be computed 
in the standard way by adding the matter action (12.11) to the Einstein-Hilbert action for the 
gravitational field. Equivalently we shall put the stress-energy tensor (12.141) in the RHS of 
the Einstein field equations of general relativity. (Obviously we can use also the gravitational 
action and field equations of any favorite metric theory of gravity.) 

Finally let us exploit some constraint equations which are satisfied by the general solution 
of the equations (12. 6p and (I2.10p . First we notice that a consequence of Eq. (I2.10p . implied 
by the antisymmetry of the Riemann tensor with respect to its first pair of indices, is 

Ufl & = 0. (2.18) 

Let us prove that in fact this relation is identically satisfied, in the sense that it is implied 
by the other equations we have. We recall that f2 M = uo p — p M . Using the definition (I2.12p 
we compute first the contraction u^uj 11 , and reduce it thanks to the easy-to-check formula 

V = -F^ ± , (2.19) 
m 

which follows from the fact that V is a function of n±/m only. In Eq. (I2.19p we have used 
the expression of the quadri-force (I2.8P and we have defined 10 

^ = ■ (2 ' 20) 



Beware of the fact that since the motion of the particle will not be geodesic (u^ 1 ^ 0) , the latter definition 
ir^_ = D[-L^ -K v ]/d,T is in general different from _L£J Dir v /dr. Thus one is not allowed to commute the 
operations of perpendicular projection _L£J and of taking the covariant time derivative Djdr. 
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At this stage we arrive at the intermediate result u^uj^ = —77^(7^ + 2F M )/(2m). Next we 
have recourse to the other equation (I2.6p . and we remind the fact that p^ = since the 
norm of p M is a constant, p 2 = Am 2 c 2 . This yields immediately u^uj^ = which is indeed 
the constraint (12.1 8p we wanted to prove. So this constraint is consistent with our basic 
equations (J2S]) and (12~TU|) . 

We now show that there is another constraint relation, which is a non-trivial consequence 
of our equations, and will be used below to find an interesting particular solution of those 
equations. This relation is obtained by contracting Eq. (12.61) with the four velocity. Because 

is orthogonal to the four- velocity [cf Eq. (12.91) ]. we obtain 

V^ = °> ( 2 - 21 ) 
which can easily be transformed, using the definition (12. 7p of P M , into 11 

An/ + (A-l)«/ = 0. (2.22) 

The constraint relation (I2.22p could be viewed as a differential equation for the quantity A 
which is defined by (12.5b[) . Since it is a consequence of our main equations (12.61) and (12.101) . 



the relation f!2.22j) is to be satisfied by any solutions of those equations. 



C. Particular solution of the equations 

The general dynamics of dipolar particles in the present approach is given by Eqs. (12. 6p 
and (12.101) . whose consequence as we have just seen is Eq. (I2.22p . From now on we shall re- 
strict our attention to a particular class of solutions owning some features that are consistent, 
as we shall see, with the intuitive idea of a dipole moment in "equilibrium" in the gravita- 
tional field. This class of solutions is obtained by solving the constraint equation f 1 2 . 2 2 j) in 
the simplest way that 

A = l. (2.23) 

This choice will greatly simplify the equations of motion (12.61) and (I2.10p . and make them 
quite attractive. We shall however leave open the question of how generic the solutions 
satisfying (12.231) are, and under which conditions (if any) could they be made unique. 12 In 
more details the relation (12.231) reads 

- u = 2mc 2 . (2.24) 
Ann V ' 



11 We can write also the alternative equivalent expression 

A(2m + u^) + A(A - l)u^ = . 

12 One could be tempted to replace A = 1 back into the original action (|2 . 1 1) , therefore defining the alternative 
action 

4m 

However this action would not describe a dipolar particle in the sense we want. Notably since as we see 
S involves now a mass term in the ordinary sense (with mass 2m), its stress-energy tensor will contain 
some unwanted "monopolar" mass contribution ~ mnu^u" . 
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This can be regarded as an equation giving the time-like component of 71"", which is parallel to 
the velocity w M (i.e. u^tc^), in terms of the space-like components, perpendicular to The 
solution of (12.241) is best expressed in terms of the vector p M (which is now = 2m w M — 7r M 
because A = 1), and we get 

u/ /c = e y/T^irHr^ Am 2 c 2 , (2.25) 

where e = ±1 tells whether p M is future or past directed. 

Using Eq. (12.231) it will become clear that the linear momentum P M represents the flow 
of inertial or equivalently passive-type gravitational mass, while the other linear momentum 
ft" is associated with some active-type gravitational mass. Furthermore, for the class of so- 
lutions satisfying (12.231) . it will happen most remarkably that the "physical" dipole moment, 
namely the one which appears in the final equations of motion and stress-energy tensor, is 
the projection orthogonal 7r" to the four velocity, rather than 7r" itself. Interestingly, we 
shall find that the longitudinal component, u M vr M , which never appears in the final equations 
and is therefore unphysical (i.e. unobservable), is actually given by a complex number, see 
Eq. (POOD below. 

From Eqs. (I2.5P and (12. 7p we find that when A = 1 the linear momentum P" is simplified 
to the following vector, which is time-like, 

F' , = 2m/. (2.26) 

Then (12.211) is satisfied simply because the quadri-norm of is constant. Thus the equation 
of motion (12.61) gives the particle's quadri-acceleration a M = w M as 

2ma" = -2P". (2.27) 

Since the mass coefficient 2m is in factor of the acceleration (which incorporates both inertial 
and gravitational effects), it can equivalently be interpreted as the particle's inertial mass 
and passive gravitational mass, M ; = M p = 2m; in this sense the equivalence principle is 
satisfied. However the motion is not geodesic as the result of the force P" which is supposed 
to reflect the internal structure of the dipolar particle. 

Let us reduce next the expression of ft", defined by (12.1 ip and (12.121) . To this end we 
make use of Eq. (I2.24p to obtain first the alternative expression 

ftM = ^ [V + u v r] + tt" - . (2.28) 

c mc 

The point is that we can express ft" entirely in terms of the orthogonal projection 7r^ =-L^ n u . 
We replace 7r" in (I2.28P by its equivalent expression in terms of the time derivative of 7r^, 
namely tc^_ = D-n^jdr already defined in (12.201) . and which we recall is different from the 
alternative object (7r M )j_ =_L^ 7r M . An easy computation, using also Eq. (12.271) . brings then 
f2 M into the simple form 

ft" = ^^+±£7^, (2.29) 
cr 

displaying the longitudinal vs perpendicular decomposition of ft M with respect to the four 
velocity. Notice that in the second term of (I2.29P the orthogonal projector appears two 
times: one explicitly in front of the term and one contained into 7r^. The expression (I2.29P 
makes it clear that 

+ V = . (2.30) 
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Motivated by the fact that fi M parametrizes the "mass term" or "monopole part" of the 
stress-energy tensor (12.141) . we define the particle's active gravitational mass M a — namely 
the mass which "actively" generates the gravitational field — as the coefficient of the velocity 

in Eq. (j2.29p . More precisely, M a is defined by the longitudinal part of f2 M along as 13 



V 



(2.31) 



This mass is not conserved because of the work done by the internal force F^, and we find 
[consistently with Eq. (12.191) 

M a = At F.W . (2.32) 



mcr 



With this notation the equation (12.291) can be viewed as the classic relation between the par- 
ticle's linear momentum f2 M and the four velocity u^, namely fi M = M a u^-\- (dipolar effects). 
A similar relation holds in the case of spinning particles, where the linear momentum and 



velocity differ by spin effects [21|, |22j, [23 



The evolution equation (I2.10p reads now 



ii 



dr 



-v n ± 



(2.33) 



where we have taken advantage of the symmetries of the Riemann tensor to replace tt u by 
7r^ in the RHS. This point is not without interest because we discover that the equation 
depends only on the orthogonal projection tc^_ — indeed recall that V itself is a function of 
tt± = Vy/^i^T- m Sect. IIII Al we shall interpret the non relativistic limit of Eq. (I2.33P as 
an equilibrium condition for the dipole moment in a gravitational field. 
Finally the stress-energy tensor for the class of solutions satisfying A 



1 is given by 



= n 



c 



c 2 v p rp^ 



where the polarization tensor (12.171) reads now 



n 

72 



u 



(2.34) 



(2.35) 



Here also we have been able to replace 7r M by tt^_, thanks to the antisymmetry of the two 
terms in (12. 35ft . 

To conclude, the choice of solution (12.231) enables one to appreciably simplify the equa- 
tions and to ease their interpretation. A consequence is that the number of independent 
components of the dipole moment is reduced from four down to three. The component of 
the dipole moment that is along the four velocity (namely u^rr^) never appears in the final 
equations and is unobservable. The physical dipole moment is entirely described by the 
orthogonal projection variable 7r^, which is a space-like vector. 



Such definition is proposed here for heuristic discussion, but is not used directly in the formalism. 
In more details this equation can also be written as [making use of Eqs. (|2.19[) and (|2.27p ] 

(V + u„nt) J + J_£ *i = -ttI R% va uPu° . 
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III. DIPOLAR PARTICLES IN THE NON RELATIVISTIC LIMIT 



A. Quasi Newtonian equations 

We investigate the non-relativistic (NR) approximation of the dynamics of dipolar par- 
ticles, as described by the solution found in Sect. Ill C[ characterized by the equations of 
motion (12.271) and (12.331) . and by the stress-energy tensor (12.34p - (12.35j) . To proceed, we 
write explicitly all factors of c's, and consider the formal limit when c — ► oo, which is equiv- 
alent to the usual v/c — > 0, where v is the typical value of the coordinate velocity of the 
dipolar particles. In the following we systematically indicate the neglected remainder terms 
and write them as some 0(c~ n ). 

We suppose that the dipolar particles evolve in the gravitational field of some ordinary 
matter system with spatially compact support. We introduce a Cartesian coordinate grid 
{£, x 1 } (with t = x°/c), and we choose it to be inertial, i.e. without rotation nor acceleration 
with respect to some averaged cosmological matter distribution at large distances from the 
local matter system. In these coordinates the metric is asymptotically flat (the local matter 
system is freely falling in the cosmological background field). 

In the NR limit the gravitational field is described by a single Newtonian-like potential 
U, whose source will be the sum of the Newtonian densities of the ordinary matter and of 
the dipolar particles. Such potential will satisfy a Poisson equation coming from the NR 
limit of the Einstein field equations. In the usual notation U has the dimension of a velocity 
squared and is of order one when c — > oo, which we denote by U = O(c ). The metric 
coefficients are given by 



This metric is accurate enough to obtain the Poisson equation satisfied by U in the NR 
approximation, and to discuss the motion of massive particles (ordinary stars) as well as 
relativistic particles (ordinary photons) in the gravitational field. 

The four-velocity of the dipolar particle is written as = (u°,u°v l /c), where v i = 
dx 1 jdt is the coordinate velocity and u° = cdt/dr. Since the particle will be non relativistic, 
v % is of order unity when c — > oo, i.e. v l = 0(c°). In the NR limit the proper time reduces 
to the coordinate time, dr = dt + 0(c~ 2 ), hence u° = c + (9(c _1 ) and u % = v 1 + 0(c~ 2 ). 
The particle's quadri-acceleration, given by a M = du^/dr + T^ a u p u a , is then found to reduce 
in the NR limit to a = O^ 1 ) and o* = d 2 x i /dt 2 + c 2 r^ + C(c~ 2 ), where d 2 x i /dt 2 is 
the ordinary coordinate acceleration. The Christoffel symbol reads Tq = —c~ 2 g l + (9(c~ 4 ), 
where g l = dJJ is the Newtonian gravitational field, so we have 



[Of course, geodesic motion simply means that a 1 = hence d 2 x i /dt 2 = g { + C(c" 2 ).] 




(3.1b) 



(3.1a) 



(3.1c) 




(3.2) 
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Consider the internal force F^, which is defined by Eq. ( 12.81) and appears in the RHS of 
the law of motion (12.271) . As a basic (and quite natural) hypothesis, we impose that this 
force exists in the NR limit, in the sense that its spatial components F % admit a non-zero 
and finite limit when c — > oo: 

F l = O (c°) . (3.3) 

Since F^ is orthogonal to the velocity, Eq. (12.91) . its zero-th component F° vanishes in the 
NR limit, F° = (9(c _1 ) [and more precisely we have F° = v t F l /c + 0(c~ 3 )]. It is now clear 
that the dipolar particle's law of motion (I2.27P reduces in the NR approximation to the 
non-geodesic equation 15 



m— = mg l -F i + 0[-) . (3.4) 



We emphasize that g l = djU represents here the total gravitational field, which is sourced 
not only by the ordinary matter but also by the dipolar particles themselves. 

Consistently with the order of magnitude (13.31) we assume that the potential function V 
from which derives the force F % is also of the same order, 

V = O (c°) . (3.5) 

Equation (13. 5p means that we preclude any constant term in the expression of V, which would 
be of the order of 0(c 2 ) and therefore be of the form Ai c 2 , where M. is some constant mass 
parameter. Such term corresponds to a mass term in the action (12.11) . and would imply that 
the dipolar particle is endowed not only by the dipole moment 7r M (and its associated mass 
m), but also by the mass M. in the ordinary sense — dark matter particles would carry 
some mass in the ordinary sense. This Ai yields an unwanted monopolar contribution to 
the stress-energy tensor (non vanishing in the NR limit), so we simply pose M. = 0. Related 
to this we note that Eq. (13.51) implies that the active gravitational mass M a = V/c 2 we 
defined in (" 12. 3 1 [) vanishes in the NR limit: 

M a = O (i) . (3.6) 

In the NR limit the dipolar particle has zero active-type gravitational mass, hence its stress- 
energy tensor is purely dipolar, in agreement with the elementary intuition from paper I. 

Consider next the other equation (12.331) . Our basic dipole moment variable is assumed 
to be such that [consistently with (I3.3P and (13.51) ] 



O (c°) , (3.7) 



hence 7r^ = 0[c *) which follows from the orthogonality to the velocity. The fact that 
■n\ vanishes in the NR limit is crucial. It implies that the four norm of 7r^ reduces to the 
Euclidean norm of tt\_ in the NR limit: n± = [dij^n^] 1 / 2 + 0(c~ 2 ). Hence the force F l 
is a function only of 7r] in this limit. 16 We need the covariant proper time derivative of 



15 Since both a and F° are of order C(c _1 ) we see that the zero-th component of the law of motion is also 
satisfied in the NR limit. 

16 More precisely we should write that F l [ir^] = F 1 ^-^] + 0(c~ 2 ), but in the following we shall identify, by 
a slight abuse of notation, the functional F l with the original one F l . 
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tt^_ which is given by tt^_ = dn^/dr + r$ p u v TTj_. Because 7r^ vanishes in the NR limit, we 
readily find that 7r^ is not coupled to gravity in this limit, in the sense that the Christoffel 
symbols make a contribution which is of higher order. Hence we find that 7r° = (9(c _1 ), and 
most importantly that the spatial components 7r]_ reduce to the ordinary time derivative, 
namely n]_ = dn\/dt + 0(c~ 2 ). The linear momentum (12.291) then becomes Q° = (9(c _1 ) 
and Q l = dn]_/dt+ 0(c~~ 2 ). Applying the same reasoning we find that Q° = (9(c _1 ) and 
Q l = d 2 7r]_/dt 2 + 0(c~ 2 ). Finally the Riemann curvature tensor in the RHS of (12.331) yields 
a coupling to the tidal gravitational field dijU. Therefore we have proved that the NR limit 
of this equation is 



The law of evolution given by Eq. ( 13. 8p is interesting for our purpose because it can be 
regarded as a condition of equilibrium for the dipole moment in a gravitational field. It 
states that when the tidal gravitational field can be neglected, the components of the dipole 
moment tt]_ stay constant (or evolve linearly in time). Note that Eq. (13.81) is formally 
identical with the equation of geodesic deviation for the (space-like) "separation" vector tt^_ 
in the NR limit. However, because of other terms in Eq. (12.331) . the equation of geodesic 
deviation will not hold outside this limit. 

Thus our equations are the law of motion (13.41) and the equilibrium condition (13. 8p . Let us 
emphasize that these equations have a structure different from those of paper I, which were 
concocted by analogy with the model of a dipole moment in electrostatics, and in particular 
violate the equivalence principle. In the present relativistic description of the dipole, we are 
consistent with the equivalence principle, and as a result the law of motion H3.4[) is different 
from Eq. (12) in paper I. We also obtain directly an equilibrium condition (13.81) for the 
dipole moment, instead of the evolution equation (13) of paper I. Amazingly, we find that in 
order to reproduce the present equations (13.41) and (13.81) of the (NR limit of the) relativistic 
model, the two RHS of Eqs. (12) and (13) in paper I should exactly be interchanged. 

As we pointed out, because of the constraint relation (12.231) . the end equations depend 
only on the orthogonal projection of the dipole moment tt^_ (or n]_ in the NR limit). Still it 
is interesting to control the original, unprojected dipole moment variable 7r M which entered 
into the action (12.11) . We do it here, by considering Eq. (I2.25P in the NR limit. From the 
results 7T° = 0(c) and n l = O(c ), we find that T^tt^t!^ = C(c°) (because of the orthogonal 
projection), hence this term is negligible with respect to the mass term in Eq. (12.251) . This 
equation thus becomes = 2eimc 2 + O(c ) where i = a/— 1 is the imaginary number 
(and e = ±1). The appearance of a complex quantity is due to the space-like character 
of p^ (which satisfies p 2 = 4m 2 c 2 ). Now we further deduce p° = —2eimc + 0(c~ l ) and 
p l = 0(c°) in the NR limit. The interesting finding is that the components of 7r M are seen 
to be complex in the NR approximation for that solution, and given by 




(3.8) 




(3.9b) 



(3.9a) 



In particular the longitudinal component along the velocity is 




2mc 2 (l + si)t + O(c ) . 



(3.10) 
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However, the components of 7r M itself do not appear into the final equations of motion 
and stress-energy tensor, and are therefore unphysical; only the components of tt^_, which 
parametrize the final equations, represent the physical variables, and these are (to be con- 
sidered as) real. 

We compute the stress-energy tensor fl2.34p - fl2.35p in the NR limit. In such limit the 
number density n of dipolar particles satisfies the Eulerian continuity equation, 

d t n + d t (nv l ) = °{^j ■ (3- 11 ) 

From Eq. (13. 5p [or equivalently (13. 6p ]. together with the fact that 7r° = C(c _1 ), we deduce 
that the monopolar term in Eq. (I2.34p is zero in the NR limit so that (most satisfactorily) 
the stress-energy tensor becomes purely dipolar. The covariant derivative of the polarization 
tensor can be approximated by an ordinary derivative, so the components of T^ u are simply 
obtained as 

T 00 = -c 2 <9aT + O(c ) , (3.12a) 
T 0l = O{c), (3.12b) 
T 3 = O(c ) , (3.12c) 

in which the density of dipole moment or polarization vector is defined by 

rP = n7ri. (3.13) 

It is clear from Eqs. (13.121) that since the T^ u of dipolar particles is to be added to the 
one of the other matter fields, the only effect of the dipolar particles in the NR limit is to 
add to the Newtonian density of the ordinary matter the density of polarization 

Ppolar = -dilV . (3.14) 

The Einstein field equations reduce in the NR limit to the Poisson equation for the potential 
U defined by the metric coefficients (13. ip . The point is that the source of the Poisson 
equation is the total matter density; therefore we have proved that in the NR limit 

AU = -4nG(p + p polar ) + O (Jz) , (3-15) 

where p denotes the Newtonian density of the ordinary (monopolar) matter while p po iar is 
the dipolar matter density found in (13.141) . 

The gravitational field, given by the metric (13.11) where U is solution of (13.151) . affects the 
dynamics of any matter distribution. In the case of the dipolar particles we have already 
derived the equation of motion in Eq. (13. 4p above. For non relativistic ordinary particles 
( "ordinary stars" assimilated as point masses) we have the standard acceleration law, coming 
from the universal coupling to gravity and the geodesic equation, namely 

dt z / ordinary \C J 

In the case of the motion of relativistic ordinary particles (photons), we find that the equation 
of motion is exactly given by what general relativity predicts, i.e. 

(^)^=r i+ ^' , -«w +0 (i)' (3 - 17 » 
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where the velocity of the photon cf3 % = (dx 1 / dt) p hoton satisfies /3 J '/? J " (1 + 2U/c 2 ) = l — 2U/c 2 . 
Indeed the result (13.171) is the consequence of the metric (13. ip 17 and the expression of the 
stress-energy tensor (13.121) . However recall that g l = diU where U is solution of Eq. (13.151) . 
so the light deflection of photons is given by the usual formula (13.171) . but in which U is 
generated by the ordinary stars and also by the distribution of dipolar particles themselves. 

The latter point on light deflection constitutes an attractive feature of the present model 
of dipolar dark matter. Indeed the distribution of dark matter is felt not only by ordinary 
stars (and gas) moving around galaxies, but also by photons, as is evidenced in experiments 
probing the large scale structure using light deflection and amplification (weak lensing ex- 
periments (13]). The problem of light deflection turned out to be crucial in the construction 
of relativistic MOND theories based on extra fields besides the metric tensor of general rela- 
tivity Il5l . |16J. To account for the observed strong light bending in lensing experiments, 



one is obliged in these theories to modify the standard conformal coupling of a scalar field 
to matter by means of some vector field (preferred or dynamical) especially designed for this 
purpose The resulting tensor-scalar-vector theory is inevitably complicated. In 



the present model, by contrast, we have derived in a natural way the formula for the light 
bending (13.171) . which takes into account, via the Poisson equation (I3.15p . the effect due to 
the distribution of dipolar dark matter. 

To summarize this Section, the equations giving the dynamics of the fluid of dipolar 
particles in the NR limit are: 

(i) Their equation of motion (13. 4ft ; 

(ii) An "equilibrium" condition for the dipole moments (13. 8ft ; 

(iii) The conservation of the number of particles or continuity equation (13. lift ; 

(iv) The Poisson equation (I3.15P for the gravitational field. 

In addition we have Eqs. (13. 16[) and (13.171) for the motion of ordinary matter and photons 
respectively. 



B. Link with the modified Newtonian dynamics 

Having arrived at the Poisson equation (13.151) . whose source contains the density of 
polarization p po i ar given by Eq. (I3.14p . we apply the arguments of paper I and recover the 
MOND equation [H, @, 0] when the polarization vector n* is aligned with the gravitational 
field. We pose 

n ' = -SG 9 *' < 318 > 
and interpret Xi following paper I, as a coefficient of "gravitational susceptibility" reflecting 
the properties of the dipolar medium. We assume that x depends on the norm g = \g l \ 
of the gravitational field, in analogy with the electric susceptibility of a dielectric material 
which depends on the norm of the electric field. Next we define 

H = l + X- (3-19) 



17 We assumed that the spatial metric is the general relativistic one, with post-Newtonian parameter 
7 = 1. Obviously the calculation could be done for a general value of the parameter 7. 
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Without loss of generality, we can scale g by means of some constant acceleration ao, so 
that (13.191) defines a certain function fi(g/ao). The Poisson equation (13.151) becomes then 
equivalent to the MOND equation (which has the form of a modified Poisson equation 18 ) 



(3.20) 



To account with astronomical observations we know that the Newtonian law (without dark 
matter) must be valid for strong enough (though non relativistic) gravitational fields, much 
above the constant acceleration scale ao- Thus /i ~ 1 when g > ao. The fact that the 
acceleration scale is the relevant one is not trivial but was found early on by Milgrom [H, 0, 0] . 
Indeed, a remarkable fit to many observations of rotation curves of galaxies has been achieved 
by assuming that p ~ g/ao in the regime of weak gravitational fields, for g <C a$. The 
acceleration scale ao is empirically found to be of the order of 10 -10 m s -2 (the same numerical 
value, of course, for all galaxies). Interestingly the numerical value of ao is close to the Hubble 
scale, a ~ cH . 

Interpolating between the Newtonian and MOND regimes we have < /x < 1, hence the 
gravitational susceptibility defined by Eq. (13. 19j) must be negative: 



This fact was interpreted in paper I within a model consisting of a pair of sub-particles with 
positive inertial masses but opposite gravitational masses ±m (in analogy with the electric 
dipole made of two charges ±g). The negative sign in (13.2 1[) was then seen to reflect the 
fact that gravity is governed by a negative Coulomb law in the NR limit — like masses 
attract and unlike ones repel [35] . As a result the polarization masses tend to increase the 
magnitude of the gravitational field, by an effect which can be referred to as gravitational 
"anti-screening", and is opposite to the usual screening of electric charges by polarization 
charges in electrostatics. The negative sign of \ was emphasized in paper I as the main 
argument for viewing MOND as a mechanism of gravitational polarization. 

In principle, our task would be to justify the proportionality relation between the polar- 
ization IT and the gravitational field g l , Eq. (13.181) . In the model of paper I we invoked an 
equilibrium in which the distance between the sub-particles constituting the dipole remains 
constant. In the present paper we have to modify the argument because the equations of 
the relativistic model (in the NR limit) are different from those of paper I. Essentially we 
shall find that when the tidal gravitational field can be neglected there is a solution for 
which the dipole moments are at rest with respect to the local matter distribution. For this 
solution the internal force F l must compensate exactly for the gravitational force. When 
this is realized we find that the dipole moment is aligned with the gravitational field. 

Suppose that the fluid of dipole moments fills an asymptotically flat space-time generated 
by the local distribution of ordinary matter. When choosing such asymptotically flat space- 
time we implicitly assume that the isolated matter system is freely falling in the cosmological 
gravitational field generated by far-distant masses in the universe. The MOND regime will 
take place in the region far from the isolated system where gravity is weak, so we expect 
that a good approximation is to neglect the tidal gravitational field: 



We are adop ting the variant of the MOND equation which is derivable from a non-relativistic La- 



-1<X<0. 



(3.21) 



dull « . 



(3.22) 



grangian 




Henceforth we no longer indicate the neglected remainder term 0(c 2 ). 
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The discussion which follows is based on this approximation. More detailed calculations 
taking into account the tidal gravitational field will be performed in Sect. IIV1 

The equation (13. 8p tells us that when (13.221) holds tt]_ is constant or varies linearly with 
time. Discarding any linear variation in time, we assume it to be constant, 

vri ss const , (3.23) 

where the approximation sign « reminds us that this is true when the tidal gravitational 
field is neglected. Now we have seen that the force F % is a function of tc]_ in the NR limit, so 
it must also be constant, F l m const. Similarly for the coordinate acceleration of the dipolar 
particle given by Eq. (13. 4p . d 2 x l /dt 2 « const. Our basic motive is that we shall assume, 
without justification if it were not by the existence of the solution, that the coordinate 
acceleration of the dipolar particles is everywhere approximately zero: 

d 2 x % 

This is consistent with the expectation that the dipole moments stay at rest with respect 
to some averaged cosmological matter distribution on cosmological scales. In this view 
the fluid of dipolar particles appears essentially to be an immobile (static) "ether" , weakly 
influenced by the ordinary matter distribution, in agreement with the model of paper I. 
Substituting (I3.24[) into the law of motion (13.41) we obtain 

Pwmg*, (3.25) 

so in this solution the internal force F l is equal to the weight mg\ exactly like in paper I. 
Notice that Eqs. (I3.24I) - (I3.25I) actually say that the particle is accelerated in a quadri- 
dimensional sense, because its motion is non-geodesic (recall ma M = — F M ). The quadri- force 
F^ acts like a rocket to accelerate the dipolar particle and hold it at rest in the gravitational 
field. 

Because F % is proportional to 7i]_ [see Eq. (12.81) ]. the result (I3.25P implies that 7i]_ and hence 
the polarization vector II 1 = niT]_ are aligned with the gravitational field. The proportion- 
ality relation (13.181) is therefore justified and we have verified the MOND equation (13.201) . 
Following paper I we then look for the expression of F % that corresponds to the specific 
MOND regime where /i ~ g/do, and obtain it in the form of an expansion when tt± — > 0: 



1 + t n± + (nl) 
a 



(3.26) 



Such expansion can be viewed as a short-distance expansion in the separation x = ir±/m 
between the dipole's constituents. We have posed k = 4irGn which is assumed here to be 
constant and uniform. 19 The potential function V from which derives the latter force reads 
(in the approximation where k is independent of 7r]_) 



2 L 



(3.27) 



Thus Eq. (|3.26p depends on the space density n of particles. Note however that (|3.26|) represents the value 
of the force at equilibrium, when the equilibrium condition (|3 . 25(1 is satisfied. This force is computed from 
the internal force defined in general by Eq. (|2.8|) . which constitutes a definition intrinsic to the dipole, 
i.e. valid for the single dipole moment independently of n. Thus, derives from the part of the 
action (|2.ip corresponding to a single particle, without reference to n. 
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It is noteworthy that this potential takes the form of an harmonic oscillator, i.e. is quadratic 
at dominant order when tt± — > (see paper I for discussion). In the present formalism the 
MOND acceleration scale a appears to be related via Eqs. ( I3.26I) - (I3.27I) to the properties 
of the internal dipolar interaction at short distances. 

To summarize, the present discussion confirms the relation between MOND and a specific 
form of dipolar dark matter. Admittedly this relation is striking - - cf the electrostatic 
analogy for the MOND equation and the polarization density (13.141) . the correct sign of the 
susceptibility coefficient (13. 2 1 p . the harmonic form (I3.27P of the potential V. Dark matter 
could consist of a fluid of dipole moments, polarized in the gravitational field of ordinary 
galaxies in the MOND regime where g < ao (but inactive in the regime of Newtonian 
gravitational fields where g 3> ao). This fluid of dark matter would essentially be static with 
respect to the averaged cosmological matter distribution at large scales. 



IV. DIPOLAR PARTICLES IN A CENTRAL GRAVITATIONAL FIELD 

The previous Section was restricted to the case where the tidal gravitational field is 
negligible, Eq. (I3.22[) . Though probably a good approximation, this restriction hides the 
way the complete equations could be integrated in more realistic situations (possibly using 
numerical methods). Recall from Sect. IIIIA1 that the dynamics of dipolar particles in the NR 
limit is given by the equation of motion (13. 4p . the equation of evolution (I3.8p . the continuity 
equation (13.111) and the Poisson equation (13.151) . In the present Section we integrate those 
equations in the particular case where the gravitational field is generated by a point mass 
M (made of ordinary matter). 

We shall assume that (1) the dipole moments are aligned with the gravitational field, and 
(2) the fluid of dipolar particles is stationary with a purely radial flow. With those hypoth- 
esis we shall basically show that there is a one-to-one correspondence between the MOND 
function fi(g/ao) and the potential function V(iT±/m) describing the "internal physics" of 
the dipole moment. Note that by (1) we assume that the dipolar medium is polarized in the 
gravitational field but do not prove it. 20 However we shall find a consistent solution which 
describes the polarization state of this medium — this partly justifies the assumption. 

The gravitational field generated by M is central. We denote by r = \x l \ the radial 
distance to the mass, by U(r) the gravitational potential, and by g % = diU = U'{r)n l the 
gravitational field, where it = x l /r is the direction to the observer and the prime denotes the 
derivative with respect to r. We have U'(r) < and g = \g l \ = —U'{r). By the assumption 
(1) the dipole moment is aligned with g l so we pose tt]_ = —ir±ri l , where the minus sign 
reflects the fact that the moment should like g % be directed towards M, in conformity with 
the sign of Eq. (I3.2ip . Here ttj_ is the Euclidean norm of for simplicity we keep the same 
name as for the four norm (12 .4p because both agree in the NR limit. Similarly the internal 
force F l , which is parallel to the dipole moment and in the same direction, is written as 
F i = -Fn\ Thus Eq. ([22]) becomes 

dV , . 

F = m- — . (4.1) 
an i 
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This could be viewed as an indication that the model is more "effective" than fundamental (like some 



models of dielectrics where the electric polarization is parallel to the electric field 19]). 
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By the assumption (2) the velocity of the dipolar particles is radial, so we write it as 
v l = dx % jdt = n % drjdt (such hypothesis will partly be justified by the consistency of the 
solution). With these notations Eq. (13. 8p becomes 

= (4.2, 

where U" is the second derivative with respect to r, while Eq. (13.41) reads 

d 2 r dV 



dt 2 diii 



g. (4.3) 



Let us now adopt an Eulerian description of the dipolar fluid. For a stationary fluid 
there is no dependence on time, and in the central potential the norms of the velocity 
and dipole moment depend only on the radial distance: v(r) and 7Tj_(r). Their total time 
derivatives are then given by the usual Eulerian derivative as dr/dt = v, then dv/dt = vv' 
and d%±/dt = vtt' ± . Equation (14. 2p is thus reduced to 

7 7r' x + V 2 7r" = 7T± U" , (4.4) 

where j — vv' denotes the acceleration, which is itself given by Eq. H4.3[) as 

dV . . 

^=d^~ 9 - (45) 

A consequence of the latter equation is 



V ' = <{j~ U ) • (46) 

We now eliminate 7 between (14.41) and (14.51) to obtain 

V + U' tt' ± + v 2 tt" = 7T ± U" , (4.7) 

which can be transformed by further manipulations [using Eq. (14.61) ] into an equation having 
the form of a "conservation" law: 

(V-v 2 7r r ± + 7r x U'Y = 0. (4.8) 

Upon integration of this conservation law we shall obtain an arbitrary constant (independent 
of r), i.e. V — v 2 7c'j_ + tt± U' = const. However this constant can be absorbed into the 
definition of the potential V (which is anyway defined up to a constant), and we arrive at 
the simple result 

V = n ± g + v 2 n' ± . (4.9) 
Moreover we have the independent equation (14.51) which we rewrite using (14 .9p as 



dV d 
9 + 



dn 1 dr 



V - TT±g 

2tt' l 



(4.10) 



The previous relations are valid for any central potential U(r), and we want now to 
specify that U(r) is generated by the point mass M. The density of ordinary matter in the 
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RHS of the Poisson equation (13.151) is therefore p = M5(x t ), where S is the three-dimensional 
Dirac function. Using also the density of polarization fl3.13p - fl3.14p . we find that the Poisson 
equation becomes 

d i {g i -kir* L ) = -ATrGM6, (4.11) 

in which we employ the notation k = 4n Gn already used in Sect. IIIIBI In the present case 
of spherical symmetry the latter equation is equivalent to 

GM 

g-k7r ± = —, (4.12) 

where the RHS is made of the Newtonian potential for the mass M. From Eq. (13.181) we 
have kn± = ~xg, where the (negative) susceptibility coefficient x is related to the MOND 
function by \i = 1 + x- Thus, 

GM 

gn = -pr, (4-13) 

which is nothing but the MOND equation in the case of a point mass source, as originally 
postulated in Refs. (1, @, 0] . 21 There is no surprise in recovering (14.131) because we have 
already shown on general grounds that the Poisson equation with a dipolar source term 
satisfying (I3.18P is equivalent to the MOND equation (I3.20p . Equation (14.131) could also be 
directly deduced from (I3.20p . On the other hand we still have the relation between and 
g which we prefer to write in terms of fi as 

k7T ± = (l-fi)g. (4.14) 

We deal next with the equation of conservation of the number of particles (13. lip . This 
equation is easy to solve because the fluid is stationary, the velocity field is purely radial, and 
we are in spherical symmetry. It suffices to say that the flow of dipolar particles crossing the 
surface of the sphere S = 4irr 2 is constant, which means that nv S = C where C represents 
the number of particles passing through S per unit of time — a constant. Combining this 
with Eq. (gZEJ) we get 

- = C (4-15) 

where C is a constant giving the flow through S per unit of central mass M (i.e. C = C/M). 
Both C and v are negative for an inward flow directed toward the central mass, and positive 
for an outward flow. We shall find that for the same configuration for tc±, n and V the 
two solutions are possible, and can be deduced from each other by a time reversal. For 
definiteness we shall choose C > 0. 

Our equations are the first integral (I4.9p . the equation of motion fl4.10p . the gravitational 
field equation (I4.13p . the polarization relation (I4.14p and the conservation of particles equiv- 
alent to (14.151) . Thus, five equations in all for the five unknowns g, ttj_, V, v and k (where 
k is equivalent to the number density n). Note that in this counting we have considered 
that the MOND function /i is known, but we have included the potential function V in the 
list of unknowns. Indeed we shall show that once /i is specified (for instance /i = 1 — e~ 9 ^ a ° 



In the extreme MOND regime we have \i ~ g/ao and the relation (|4.13|) yields the famous logarithmic 
potential U ~ —y/GMao lnr explaining the flat rotation curves of galaxies and the Tully-Fisher relation 

^flat ~ (GMfflo) 1 / 4 . 
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as chosen below) the potential function V is determined. It would seem a priori more nat- 
ural to proceed the other way around, namely to specify first V because it is part of the 
more fundamental action (12.11) . and only then to deduce \x which would tell which kind of 
MOND phenomenology this action corresponds to. Since there will basically be a bijective 
correspondence between V and /i, it is clear that the two approaches are equivalent, and 
we can say that a choice for the action (12.1 1) determines jj, as expected (at least under the 
hypothesis of the present calculation). 
From Eqs. (I4.14I) - (I4.15I) we obtain v as 



V = C 7T I 



Replacing it into (14. 9 p yields the differential equation 



(4.16) 



7T 



C 2 n ± 



'V 

9 

7T_l 



Furthermore we can write (14.101) into the equivalent form 

V V 2 , /i// 

7l ± 7T_l (1 - jl) 



(4.17) 



3 ' 



(4.18) 



Since fi is given, g is known by Eq. (I4.13p . and the equations (14 . 1 T[) — (14 . 1 8 j) form a system of 
coupled differential equations for ir± and V. In a first stage the ratio between V and ir± is 
integrated with the result 

/i 



V 

7lj_ 



9 + 



9 



dg 



(4.19) 



(The standard notation is used for the indefinite integral, defined modulo an integration 
constant.) Then the solution for n± is found to be 



"" 1 = (^/ rfri 7 



1/2 



(4.20) 



The velocity field is then given by (14.161) . the density of dipole moments n — hj (47rG) follows 
from (14.151) . and V is deduced from (I4.19p . As we said, for any choice of function /i we can 
determine V — and vice versa. In this sense we have related the phenomenology of MOND 
to some more "basic" physics associated with the description of dark matter. 

Let us exemplify the previous resolution by showing the case of the MOND function that 
corresponds to the simple susceptibility coefficient x — —e~ 9 ^ a °, namely 



1 - e- 9/ao . 



(4.21) 



This function is obviously interesting (though maybe very special) because it is exponentially 
close to its Newtonian limit when g —>■ oo. In this case we readily integrate Eq. (14.191) and 
find 

_ = ge g/a _ ao ( e s/ao _ ^ ln ( e s/ao _ X ) _ (422) 
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DIPOLE MOMENT VARIABLE 




0.5 1 1.5 2 2.5 3 

FIG. 1: Dipole moment tt± as a function of the rescaled distance p = r/r^ denned by (|4.25p . On the 
horizontal axis p = 1 corresponds to the transition between the Newtonian and MOND regimes. 
The unit of dipole moment on the vertical axis is (GMao) 1//4 jC . 



Here the integration constant has been chosen in such a way that the ratio V '/tt± becomes 
equivalent to g in the Newtonian regime g — > oo (any other choice would make it growing 
exponentially). The dipole moment is then given by 

/ 2a f r , ,ln(l-e-^o)\i/2 

where g' = g(r'). We have fixed the integration constant (in the lower bound of the integral) 
in order to ensure that tcj_ — > in the Newtonian limit corresponding to r — > 0. Using the 
fact that g(r') is implicitly given by Eq. (14.131) . we obtain 

(CM*)"' ( ^__ 3/2 __„! + (» -l)e~",_„ __„^ 1/2 



In Fig. [1] we plot the dipole moment 7r^ as a function of the rescaled distance defined by 
p = r/r Q , where r is the characteristic length at which the transition between the Newtonian 
and MOND regimes occurs, i.e. which is such that g(ro) = ao and is given by 



GM 



1/2 



\a (l-e v )j 

In Figs. [2] and [3] we plot the velocity v and the number density n of dipolar particles as 
functions of p = r/rQ. We give the polarization of the medium in Fig. HI the density of 
polarization masses in Fig. [5l and the potential function V in Fig. [61 

As we can see from Fig. [T]the moment 7r^ tends (exponentially) to zero when r — > 0, so 
the dipole moments actually do not exist in the vicinity of the ordinary mass. However 7r^ 
diverges in the MOND domain where r — > oo; more precisely, we can check that ttj_ behaves 
like ~ r(lnr) 1 / 2 when r — > oo. Moreover, we see from Fig. [3] that the number density n of 
dipole moments decreases at infinity, at the rate ~ r _2 (lnr) -1 / 2 . To ensure the conservation 
of the mass flow the velocity v must then increase like ~ (lnr) 1 / 2 , which is indeed what we 
observe in Fig. [2] (recall that the particles do not obey the geodesic equation). 
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VELOCITY OF DIPOLAR PARTICLES 
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FIG. 2: Velocity of the dipolar particles v as a function of the distance p. The unit of velocity on 
the vertical axis is (GMao) 1 ^ 4 . 

NUMBER DENSITY OF DIPOLAR PARTICLES 
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0.5 1 1.5 2 2.5 3 

FIG. 3: Number density of the dipolar particles n as a function of the distance p. The unit of 
number density on the vertical axis is (a^/GM^^C /(4irG). 

The fact that the dipole moment ir± opens up more and more and diverges when r 3> tq 
(deep in the MOND regime, far from the ordinary matter) may appear to be surprising. 
However, what is important is the density of dipole moments or polarization II = nir±. As 
we have seen the number density of dipole moments tends to zero in the MOND domain 
r — > oo, and as a result the polarization (plotted by Fig. 0J does decrease at large distances 
as expected (at the rate ~ r _1 ). 22 On the other hand, we see that with the present choice 
of MOND function (I4.2ip the polarization decays exponentially in the Newtonian regime at 
short distances. 

In the present model the polarization reflects the distribution of the gravitational field 
that is induced by the fluid of dipolar particles. The polarization shown in Fig. H]is in fact 



Figure [4] is the analogue of Fig. 1 in the Newtonian- like model of paper I (indeed the density n of dipole 
moments is considered to be constant in paper I, so the dipole moment varies like the polarization). 
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POLARIZATION OF THE MEDIUM 

0.4 r 



. 3 



0.2 



0.1 




2 4 6 8 10 



FIG. 4: Polarization of the medium II as a function of the distance p. The unit of polarization on 
the vertical axis is ao / (47rG) . 



DENSITY OF POLARIZATION MASSES 




FIG. 5: Density of polarization masses p po iar as a function of the distance p. The unit of mass 
density on the vertical axis is ao/(inGro). 



given by the analytic formula 

(Notice that this formula, equivalent to (I3.18p . could be written directly from the standard 
MOND equation [1, @, 0] ; only the interpretation we propose here as the polarization of the 
dipolar medium is new.) 

Even more relevant than the polarization to illustrate the effect of dark matter is the 
density of the polarization masses p po iar defined by Eq. (13.141) . In the present case [cf the 
assumptions (1) and (2) above] we find that 

2n 

Ppoiar = II' + — . (4.27) 

r 
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POTENTIAL FUNCTION V 




12 3 4 5 



FIG. 6: Potential function V as a function of the distance p. The unit of V on the vertical axis is 
(GMalY^/C. 



Figure [5] displays p po iar as a function of the distance p = r/r . In the present model, p po iar 
really represents the mass density of dark matter (namely, that density which is added to the 
density of ordinary matter in the source of the Poisson equation for the Newtonian potential 
U). We can check from Fig. [5] that p po iar has the correct behaviour ~ r~ 2 at infinity which 
reproduces the flat rotation curves of galaxies; in more details we have 



47rr 2 V G 
which corresponds to the polarization mass 



Ppoiar ~ ~7~~~~o \l ~pr~ when r -> oo , (4.28) 



M po i ar ~ J r when r -> oo . (4.29) 

Note also that a prominent feature of Fig. [5] is that it predicts no accumulation of dark 
matter in the close vicinity of ordinary masses, when r <C Tq. Of course we meet here the 
natural explanation by MOND (and therefore also by the present model) for the absence of 
observed cusps of dark matter in the central regions of galaxies. 

We have thus shown how the equations of motion and evolution of the dipolar particles 
can be integrated in a specific example. Although the solution we have considered is rather 
idealized — the fluid of dipolar particles is stationary and the dipole moments are exactly 
aligned in the central gravitational field of a point mass 23 — it has the merit of exhibiting in 
details the link between the phenomenology of MOND and the physics of dipolar particles, 
specified by the potential function V(ir±/m) which enters the action 02. ip [the function V 
is plotted in Fig. [6] for the particular case of MOND function (I4.21I) ]. 



23 Note that the question of the stability of this solution against gravitational perturbations has not been 
investigated. 
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V. SUMMARY AND CONCLUSIONS 



Motivated by the quasi- Newtonian model of paper I [181 ]. we interpret the phenomenology 
of MOND as resulting from an effect of gravitational polarization of a medium made of 
dipole moments aligned in the gravitational field of ordinary masses. We propose an action 
principle, based on the matter action (12.11) . to describe the dynamics of dipolar particles 
and the evolution of the dipole moments in standard general relativity. The action involves 
a kinetic-like term for the evolution of the dipole moment and a scalar function V which is 
supposed to describe (at some effective level) the non-gravitational interaction between the 
constituents of the dipole. 

The dynamical variables are the dipolar particle's space-time position and the dipole 
moment four vector 7r^, which are varied independently yielding the two basic equations (12. 6p 
and (I2.10p . Variation with respect to the metric yields the stress-energy tensor (12.141) . We 
find a particular class of solutions, defined by the constraint (I2.23p . corresponding to the 
intuitive idea of a dipole moment in equilibrium. For this class of solutions the physical 
dipole moment variable is the projection tt^_ orthogonal to the particle's four velocity. 

The non-relativistic (NR) limit of the model is investigated next. The basic equations 
in this limit are (i) the equation of motion H3.4[) of the dipolar particle, (ii) an equilibrium 
condition (I3.8P for the dipole moment 7r^, (iii) the conservation law (13. lip of the number 
of particles n, and (iv) the Poisson equation (13.151) for the gravitational potential U . The 
equations (i) and (ii) are different from those of the quasi-Newtonian model of paper I 
(indeed the model of paper I violates the equivalence principle and cannot result from the 
NR limit of a general relativistic model). In addition we have the equations of motion of 
ordinary massive particles (13.161) and of photons (I3.17p . 

Assuming that the tidal gravitational field dijU can be neglected, we find that there is 
a solution for which the dipole moments are aligned with the gravitational field. For that 
solution the Poisson equation for the gravitational field reduces to the MOND equation like 
in paper I, confirming the close relation between the phenomenology of MOND and the 
dipolar dark matter. 

Finally the equations in the NR limit have been integrated in the idealized case where 
the fluid of dipolar particles is stationary and the dipole moments are exactly aligned in the 
gravitational field generated by a point mass M. The polarization of the medium is plotted 
in Fig. H] and the density of polarization masses (which represents in this model the density of 
dark matter accumulated around M) is shown in Fig. [5j There is a correspondence between 
the potential function V in the action and the MOND function fi. 

To conclude, we usually face two alternatives to the issue of dark matter: Either accept 
the existence of cold dark matter particles but which fail to reproduce in a natural way the 
rotation curves of galaxies, or postulate an ad hoc alteration of the fundamental theory of 
gravity (MOND and its relativistic extensions). In the present paper (following paper I) 
we proposed a third alternative: Keep the standard law of gravity but add to the ordinary 
matter some non-standard dark matter in order to "explain" MOND. More precisely we 
invoke a mechanism of gravitational polarization, in which the ordinary masses (galaxies) 
are "anti-screened" by polarization masses associated with gravitational dipoles playing the 
role of dark matter. The dipolar dark matter particles are in equilibrium in the gravitational 
field because of their internal structure driven by some postulated non-gravitational force 
(whose fundamental origin is unknown). 
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APPENDIX A: VARIATION OF THE ACTION FUNCTIONAL 

In this Appendix we show how to vary the dipolar action functional (12.11) . which is of the 
general form 



dr L 



(Al) 



The dynamical variables are the particle's space-time position x p (t) and the dipole moment 
carried by the particle tt p (t), both depending on the proper time parametrizing the world 
line and denoted by dr = —g pv dx p dx v / c 2 . The Lagrangian L is a function of the dy- 
namical variables through the metric g pv {x) evaluated at x p (t), the particle's four velocity 
u p = dx p /dr, the dipole moment itself vr M (r), and the covariant time derivative of the mo- 
ment 7r M = dir p " /dr + Y p u p/ n a . The four variables g^ u , u p , ir p and 7r p are considered to be 
independent in Eq. (j All) . 

We first establish a relation stating that the action ([AID is a scalar vis-a-vis the arbitrary 
infinitesimal coordinate transformation x' p = x p + e p (x). The requested linear transforma- 
tions of the independent variables are 



u 



'I'- 



ll! 



u p d n e p 



>v = * n p + -K p d e p 



7r' p 



} + 0(e 2 ) 



n p + n p d p e p 



(A2) 



together with dr' = dr which is a scalar. Here the prime refers to the new values of the 
components of the tensor when evaluated at the same space-time event. In the last line we 
have used the fact that tt p is a four vector. The "scalarity" of the action means that, under 
the previous coordinate transformation, 



dr L 



dr' L 



2*7, 



dL 



MP 



dg, 



u 



up 



dL 
du p 



dL 

d-K p 



+ 7T" 



dL 

dn p 



(A3) 



in which the same functional L appears on both sides of the equation. Using Eqs. (1A2I) we 
immediately deduce the scalarity condition 



(A4) 



We first vary the action with respect to the moment 7t p (t) subject as usual to the condition 
that tt p (±oo) = 0. The variations are 5ti p and 5tt p = d5n p / dr J r T p a u p 8it eT . We integrate the 
ordinary time derivative dSw^/dr by part, and find that the role of the Christoffel symbol 
is to make the corresponding field equation covariant, 



D f dL 



dr \ d-hv 



dL 

d-K^ 



(A5) 
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where D/dr is the covariant time derivative. The variation with respect to the position x M (r) 
[satisfying x M (±oo) = 0] is more involved, due to the dependence of dr, and g pv on the 
position. Many Christoffel symbols and their derivatives are generated in the calculation. 
The partial derivative of the Lagrangian with respect to the metric is simplified with the 
help of the condition (1A4D . We finally obtain a manifestly covariant equation, given by 



D_ f dL 

dr \d^ + Ul1 



u 



duP d-kP 



L 



dL 

dk u 



(A6) 



where R v n ,,„ is the Riemann curvature tensor. 



Finally we obtain the stress-energy tensor T pv of the fluid of dipolar particles with number 
density n satisfying the continuity equation 



(ra/ 



0. 



(A7) 



We thus vary the action with respect to the metric, which enters explicitly in the first slot of 
L in Eq. (1A1I) . and implicitly through the proper time dr and the covariant time derivative 



t\P of the dipole moment 



21 



n 



In this way we obtain 
.dL 



2 £L + 



dg r 



+ V, 



n 



u'- 



duP 



+ 7T 



dL 

d-k„ 



dL 



91 L 






) 


dkP 






dL 


) 


dk v 



(A8) 



in which we denote dL/dk p = g pu dL/ dk a . The second, dipolar- type term in this expression 
takes the form of the divergence of a tensor, and has a structure similar to some related term 
in the Belinfante-Rosenfeld 36|, |37|] symmetric stress-energy tensor of integer spin fields. An 
alternative expression of (1A8j) is readily derived with the help of the equation of motion (1A5|) 
and the scalarity condition (1A4I) ; we have [with dL/du p = g pa dL/du cr ] 



n u 



V, 



dL 

dUy) 



u 



n 



dL 



7T 



dk, 



— 7T 



dL 
duP 
dL 



+ 7T' 



wP 



dL 

dkP 



L 



(/-> 



dk 



PA 



u 



We check that the stress-energy tensor is conserved "on shell", i.e. 
motion ( 1A5I) and ( 1A6I) are satisfied: 

V„T^ = . 



(A9) 

when the equations of 
(A10) 



Finally let us make the link with the linear momenta P M and introduced in Eqs. (12. 7p 
and (12. lip , and which we have seen characterize entirely the dipolar particle. These momenta 
appear in fact to be the conjugate momenta associated with the dynamical variables n p and 
x p respectively, in the sense that 



2m 



dL 

dkP 



(Alia) 



24 The variation of the Christoffel symbol present in ttP is a tensor given by the Palatini formula as ST^ 
\gP v [V p 5g av + V a 5g pv - V v 5g pa ]. 
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it' 



9L 



+ 7T 



9^ 



(Allb) 



The relation for f2 M seems to be more complicated than a simple conjugation relation, but 
this is due to our use of the parametrization by the proper time dr in the action, instead of 
a parametrization which is independent of the dynamical variable x^. With such definitions 
the equations of motion flA50 and (IA60 become 



P„ 



2m- 



dL 



2m 



(Al2a) 
(Al2b) 



and are seen to be equivalent to the expressions derived in Eqs. (12.61) and f)2.10p [indeed 
F M = —mdL/dirt 1 with the explicit action (12.11) ]. Furthermore, using the alternative expres- 
sion (1A9|) . the stress-energy tensor is obtained as 



T" u =nQ^u" ) V, 

2m f 



n 



71 



U 



(A13) 



which is precisely the result given by Eq. (I2.14p . 
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